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Abstract. - We illustrate the Jarzynski equality on the exactly solvable model of a one- 
dimensional ideal gas in uniform expansion or compression. The analytical results for the 
probability density P(W) of the work W performed by the gas are compared with the results 
of molecular dynamics simulations for a two-dimensional dilute gas of hard spheres. 



Exactly solvable models play an important role in statistical mechanics. They complement 
and verify results that are derived from general and usually abstract arguments, while at the 
same time they offer insight and intuition. They are particularly useful in far-from-equilibrium 
situations, for which few generic exact results are available. In this letter, we focus on a 
remarkable result in nonequilibrium statistical mechanics, namely the Jarzynski equality [1,2], 
which has given rise to a certain amount of confusion about its validity and interpretation [3] . 
The Jarzynski equality relates the statistics of the amount of work W performed by a system 
in a nonequilibrium transition between two equilibrium states, to the difference in the free 
energies of these states. More precisely, one has: 



where W is the work delivered by the system upon varying an external control parameter 
following a specified schedule between an initial and final value, starting from an initial state 
of the system sampled from a canonical distribution at temperature T (/3 = fc^T). W 
is a random variable due to the sampling of the initial state. The first surprise is that the 
above specified average (...) with respect to W is independent of the schedule according to 
which the control parameter is changed between the specified initial and final values. In 
particular it is independent of whether this schedule keeps the system close to equilibrium 
(quasi-static transformation) or whether it entails large deviations from equilibrium. As a 
consequence the average is expressed in terms of the difference AF in free energy of the 
canonical equilibrium states at temperature T at the final and initial values of the control 
parameter, respectively. The second surprise is then that the equilibrium quantity AF can be 
obtained by an ensemble average over nonequilibrium measurements. Finally and foremost, 
the appearance of an equality in far from equilibrium dynamics is very surprising. In fact, the 



(exp(f3W)) = exp(-(3AF) 



(1) 
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equality Q leads, upon application of the Jensen inequality, to the familiar inequality, (W) ^ 
W rev = —AF, corresponding to the formulation of the second principle of thermodynamics 
for a system in contact with a heat bath. 

Our purpose here is to verify and complement the discussion of the Jarzynski equation 
by deriving the analytic expression of the probability density P(W) for the work W in a 
system with Newtonian dynamics. We will consider the so-called Jepsen gas, for which a 
number of other exact equilibrium and noncquilibrium results have been derived [4]. The 
Jepsen gas consists of N identical point particles of mass m moving on a line and undergoing 
perfectly elastic collisions. Actually, since the speeds are merely exchanged upon collision, 
and the identity of the particles is irrelevant, the same model can represent an ideal gas 
in which particles do not interact with each other. As initial state, t = 0, we consider a 
system at thermal equilibrium at temperature T, i.e., the particles are uniformly distributed 
in the interval [— L,0], and their velocities are randomly and independently chosen from the 
Maxwellian distribution tp(u) — \J m[3/2iT exp(—(3mu 2 /2). The right hand side of the interval 
is formed by a piston of infinite mass, which is moved according to a specified schedule. 
Although more complicated situations can be considered, we will concentrate on the case 
of a piston moving at a constant velocity V (positive or negative, i.e., corresponding to gas 
expansion or compression, respectively), from the initial position X = to the final position 
X = Vt. A similar model has been discussed in [5]. Whereas the latter paper focuses on the 
case where the work distribution is dominated by correlated recollisions of a single particle 
with the piston, we will study the so-called thermodynamic limit of an infinite system, with 
L — ► oo, N — > oo and fixed density n — N/L. In this limit, recollisions of particles with the 
piston are no longer possible (at least not for the schedule considered here). Furthermore, as 
we will show below, a comparison with the results for real but dilute gases becomes relevant. 

To evaluate W, we first note that a particle of velocity Uj > V colliding with the piston 
will recoil with the velocity — 2V — it,. Hence there is an energy transfer in the interval 
[0,t] from this particle to the piston given by AWi = 2mV(ui — V)6{xi + uit — Vt), where 
the Heaviside function expresses that the collision between the particle (with initial position 
Xi < and velocity m) and the piston has to take place before time t. We conclude that 
the probability to have a total energy transfer W from the gas to the piston during the time 
interval [0, t] is given by: 



N 



N I r°° dk { 

P(W) = (S(W-J2 A Wi)) = ^J_ ^expj ? 

/OO 1 T 

exp(ikW) (expi-ikAWj)}" 
-oo ^ 7r 

= J — exp^ikW-nt J du (u - V)<j){u) [l - exp {~2ikmV{u - V)) } j 



(2) 



The average (. . .)o is taken over the distribution of the initial positions and velocities of the 
particles. Note that we have not taken into account recollisions of the particles with the 
piston, hence only the final line, in which the thermodynamic limit has been taken, gives the 
exact result for P(W). 

The schedule under consideration can be conveniently characterized by the following two 
dimcnsionless variables: 



v = V(f3m/2) 1/2 , t = nt/{2m(3) 1/2 , 



(3) 
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Fig. 1 - Profiles of the nonsingular part of P(w) for various velocities v of the piston at a fixed 
time t = 5 (left: gas compression; right: gas expansion). The area below the curves represents the 
probability to have a nonzero value of w and is equal to [1 — exp(— tCq)]. The graphs illustrate the 
agreement between analytical calculations (solid lines) and molecular dynamics simulations results 
(dashed lines), see the main text. 



which are, essentially, the velocity of the piston measured in terms of the thermal speed of 
the gas particles, and the average number of collisions during the considered time interval 
[0, t] for a stationary piston. In terms of the scaled work w = f3W , the expression for the 
probability density becomes: 

P(w) = J°° ^exp[iqw-rC(q)} , (4) 

where the function C(q) can be written as the sum of two parts, C(q) = Co + C(q). The 
"collisionless" part, 

Co = l/\/7r exp (— v 2 ) — werfc(w) , (5) 

(where erfc(...) is the complementary error function) corresponds to the absence of collisions 
between the gas particles and the piston, i.e., when no work is performed. This leads to 
a singular contribution exp(— tCq) S(w ) in the expression of P(w). This contribution is 
exponentially decaying with respect to the scaled time r, while the damping exponent Cq is 
a rapidly decaying function of v. The second part, 

C{q) = u(l + 2iq) erfc [v(l + 2iq)} exp [v 2 {\ + 2iq) 2 - v 2 ] - 1/Vtt exp (-v 2 ) (6) 

is determined by the collisions between the gas particles and the piston, and leads to a mono- 
modal nonsingular contribution to P(w), as illustrated in Figs. ^ 121 and 01 

Based on this exact result, we first turn to the explicit verification of the Jarzynski equality. 
One has: 

f°° dq f°° 
(exp(,£W)) = (exp(iu)) = / — exp(— rC(q)) / dw exp(iqw + w) 

J — 00 ^ J —00 

= exp(-rC(i)) = exp(2rw) = exp(ntV) . (7) 
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Fig. 2 - Profiles of the nonsingular part of P(w) at different times r for fixed velocities of the piston, 
v — 1 (left) and v = —0.1 (right). Note the evolution from a highly asymmetric profile at short times, 
Eq. 114L to an asymptotic Gaussian shape with mean and standard deviation given by Eqs. 110H and 
lllll respectively. 



On the other hand, the change in the equilibrium free energy of the ideal gas due to the 
variation of its volume is a purely entropic factor given by: 



exp(-/3AF) 



L + Vt 
L 



N 



thermodynamic limit exp(ntV) 



(8) 



so that the Jarzynksi equality (yi is indeed reproduced. 

We next examine the characteristic properties of P(w). While an explicit evaluation of 
the Fourier transform appears to be difficult, one can easily obtain exact results for the 
moments using the characteristic function of the probability density, namely 



G(q) = (exp(-iqw)) = 



(— ia) n 



(9) 



In particular, one obtains the mean value of the transferred energy 

(to) = 2vt [(1 + 2v 2 )erfc(v) - 2/^/ttv cxp (-v 2 )] , 



(10) 



that increases linearly with time and has a maximum as a function of v (see Fig. |3J| . The 
centered moments of second through fourth order are 



cr 2 = (w 2 ) - (w) 2 = 8v 2 t 



(1 + v 2 ) exp (-v 2 ) - v(3 + 2u 2 )erfc (v) 



({w ~ (w)) 3 ) = 16v 3 



(3 + 12v 2 + 4v 4 )erfc(w) - (5 + 2v 2 ) exp (-v 2 ) 



(11) 



(12) 



(O- (w)) 4 ) =3f7 4 + 64 V 4 r 



2/V^F (4 + 9v 2 + 2v 4 ) exp (-u 2 ) 



- v (15 + 20v 2 + 4v 4 )erk(v) 



(13) 
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Fig. 3 - The mean value of the energy transfer per unit time {w)/t, and, in inset, ({w) — w rev )/r, 
both as a function of v. Note that (w) ^ w rev , in agreement with the second law of thermodynamics. 

Table 1. — Comparison of the mean values and higher moments obtained from the molecular dynamics 
simulations with analytical results (values in parentheses), for different values of the velocity v of the 
piston and a fixed value of |u> re „| = 2\v\t = 20. The agreement is satisfactory, except for a few cases, 
cf. the discussion in the main text. 



The explicit analytical expression of P(w) can be derived in specific limits. 

(i) In the limit of large number of collisions r ^> 1, the work is the sum of a large number 
of independent contributions. One finds, in agreement with the central limit theorem, that 
the distribution function P(w) converges to a Gaussian. The mean and standard deviation 
are given by Eqs. JTHJ and (jnj. The skewness 73 = {(w - (w)) 3 )/a 3 = 0{t- 1/2 ) and the 
kurtosis 74 = {{w — (w)) 4 )/a 4 — 3 = C(t _1 ) decay to zero. 

(ii) For the limit of small number of collisions t < 1, the moments generating function 
reads G(q) p=> 1 — rC(q), hence 



P(w) pa (1 - tC )6(w) + 



\w\ exp 



/ w y 



9(wv) , 



(14) 



corresponding to an asymmetric probability density profile (compare, for v > 0, with Eq.(18) 
from Rcf. [5] for the case of a single particle, n — l/L). 

(iii) The quasi-static limit. The reversible work W rev delivered by the ideal gas during an 
expansion over a distance Ax is PAx, with the pressure P given by rife^T. This result has to be 
compared with the irreversible expansion at a finite speed V over the same distance, i.e., Ax = 
Vt. In terms of the previously introduced dimensionless variables, this implies w rev = 2vt . 
This result is indeed recovered in the limit \v\ <C 1, since C(q) ps 2viq, and P(w) ps 5(w—w rev ). 
By including a first order correction, C(q) ps 8v 2 q 2 j \fn + 2viq (l + 4v/v / 7r) , the probability 
density P(w) assumes a Gaussian profile, with mean and variance: 



8tv< 



(w)'< 



16tv' 



(15) 



Note that the Jarzynski equality Eq. Q implies that for a Gaussian distribution the 
fluctuation-dissipation ratio 

*-<■'> -':»:, (i6) 



2{w rev - (w)) 

is equal to 1. This is indeed the case in the quasi-static limit, cf. Eq. (|15f) . but not in the 
long time limit, cf. Eqs. (fTUf> and i|TTj) . The origin of this problem can be traced back to 
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Fig. 4 - Profiles of the nonsingular part of P(w) for a fixed value of |uv e „| = 2\v\t = 20 and different 
velocities v of the piston (left: compression, right: expansion). The inset shows the same plots for 
faster compression. The arrows indicate the theoretical means, Eq. 1101 . The theoretical results (solid 
lines) agree qualitatively with the computer simulation (dashed lines) . Quantitative disagreement for 
fast compression and slow expansion are presumably due to recollisions and finite size effects. 

the contribution of large deviations to the Jarzynski average: this average corresponds to 
the characteristic function evaluated at the complex unit, G(q = i), cf. Eqs. 10 and 
Even when higher order cumulants converge to zero for the rescaled variable (w — (w))/a 
(central limit theorem) , there is no guarantee that the contribution of the non-Gaussian tails 
can be neglected in the Jarzynski average. Hence the application of the Jarzynski equality 
using a Gaussian ansatz is not reliable and can lead, in a numerical or real experiment where 
a Gaussian distribution is observed, to an erroneous value of the corresponding free energy 
difference. 

We finally turn to the practical and experimental relevance of the Jarzynski equality. In 
particular one may wonder to which extent the above ideal gas results are representative of 
experimentally accessible measurements in a dilute gas. As a first step in answering this ques- 
tion, we have performed extensive molecular dynamics simulations of a dilute two-dimensional 
gas of hard disks. Note that such simulations allow to investigate parameter regions of high 
piston speed and short times which may be difficult to reach in experiment. The obtained 
simulation results represent an average over a half million runs, for a dilute hard disk gas with 
N = 2000 disks (diameter d — 1, mass m = 1). The initial positions and velocities of the 
disks are sampled from a microcanonical ensemble in a cylinder of length L = 10 4 and cross 
section S = 10 2 (i.e., initial gas density p = 0.002), and initial "temperature" T = 1. Extra 
caution was taken to reduce correlations between samples. To compare the results with the 
one-dimensional model, the effective projected density n — pS has to be used. 

A systematic comparison between analytical and simulation results was performed for a 
wide range of values of the speed v and of time t, cf. Figs, ^and El for an illustration. Overall, 
the analytical and the simulation results agree qualitatively. In particular, the progressive 
change in the general shape of the probability distribution from adiabatic to non-adiabatic 
regime is well reproduced by the numerics. The comparison of moments given in Table 1 
confirms this agreement. There are, however, some notable discrepancies, namely for the 
mean and variance for the compression case (v < 0) and the variance for the slow piston 
(\v\ = 0.01 in Table 1 and Fig. 0J. These deviations can be explained by the fact that the 
molecular dynamics simulation deviates in two basic assumptions from the analytical model, 
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namely the ideality of the gas and the thermodynamic limit. Although a rather low density 
is used, recollisions of the gas particles with the piston are not negligible. In particular, they 
are clearly visibly in the simulations when the piston compresses the gas. This non-ideality 
presumably causes the shift in the mean value and the variance for v < 0. For very slow 
piston motion, recollisions are not a major problem, but the absence of the thermodynamic 
limit causes discrepancies between theory and simulations. Since r = 1000 for |v| = 0.01, 
several hundreds to a thousand particles collide with the piston out of the limited number of 
N=2000. Clearly the absence of fast particles in the tail of the Maxwellian distribution will 
cause a narrowing of the observed P(w) distribution. 

In conclusion, the verification or exploitation of the Jarzynski equality to measure free 
energy differences appears to be quite intricate. Due to the role of extreme events, the gain 
in computational time, when using a fast route between initial and final values of the control 
parameter, is outdone by the exponential increase in the required statistics [6]. The above 
explicit calculation and molecular dynamics simulations provide a dramatic example of this 
problem: even in the limit when a Gaussian distribution correctly describes the whole proba- 
bility mass (in the sense of the central limit theorem), extreme events may still be needed to 
correctly perform the Jarzynski average. In particular, the fluctuation-dissipation ratio 
need not approach the value 1 in a limit where P(w) converges to a Gaussian distribution. 
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